Mon. Not. R. Astron. Soc. 000, [T]-[TO] (2010) Printed 13 April 2011 (MN WI^ stylo file v2.2) 



The Intrinsic Fractions and Radio Properties of Low 
Ionization Broad Absorption Line Quasars 

Xinyu Dai^*, Francesco Shankar^, and Gregory R. Sivakoff^ 

IHomer L. Dodge Department of Physics and Astronomy, University of Oklahoma, Norman, OK 73019, USA 
7— H 1 2 Max- Planck- Institut fiir Astrophysik, Karl-Schwarzschild-Str. 1, D-85748 Garching, Germany 
^— H ' ^Department of Astronomy, University of Virginia, Charlottesville, VA 22904, USA 

o 

(N 

<C ■ ABSTRACT 

I Low-ionization (Mg li, Fe II, Fe III) broad absorption line quasars (LoBALs) probe 

a relatively obscured quasar population, and could be at an early evolutionary stage 
for quasars. We study the intrinsic fractions of LoBALs using the SDSS, 2MASS, 
and FIRST surveys. We find that the LoBAL fractions of the near infra-red (NIR) 
' and radio samples are approximately 5-7 times higher than those measured in the 

I optical sample. This suggests that the fractions measured in the NIR and radio bands 

are closer to the intrinsic fractions of the populations, and that the optical fractions 
are significantly biased due to obscuration effects, similar to high-ionization broad 
absorption line quasars (HiBALs). We also find that the LoBAL fractions decrease 
Q . with increasing radio luminosities, again, similar to HiBALs. In addition, we find 

' tentative evidence for high fractions of LoBALs at high NIR luminosities, especially 

C/3 , for FeLoBALs with a fraction of ~18 per cent at Mk^ < —31 mag. This population 

^ , ■ of NIR luminous LoBALs may be at an early evolutionary stage of quasar evolution. 

We use a two-component model of LoBALs including a pure geometric component 
. and a luminosity dependent component at high NIR luminosities, and obtain better 

^ ' fits than those from a pure geometric model. Therefore, the LoBAL population can 

I be modelled as a hybrid of both the geometric and evolutionary models, where the 

■ geometric component constitutes 3.4±0.3, 5.8±0.4, and 1.5±0.3 per cent of the quasar 

' population for BI-LoBALs, AI-LoBALs, and FeLoBALs, respectively. Considering a 

, population of obscured quasars that do not enter the SDSS survey, which could have 

T^lj- ■ a much higher LoBAL fraction, we expect that intrinsic fraction of LoBALs could be 

' even higher. 
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1 INTRODUCTION Fe ii; LoBALs). The majority of BALQSOs are HiBALs. 

/T^ . T 1 1 In fact, all LoBALs also contain the high-ionization troughs 

Broad absorption line quasars (BALQbOs) are a sub-sample . . , , , , „ , , 

^ ' m their spectra (e.g., Weymann et al. 1991; Irump et al. 

or quasars exhibiting blue-shiited absorption troughs (e.g., nnnn\ n • i ii r i ■ • ^ i 

N T 1 fin rTTr 2006). Besioes the presence oi low-ionization troughs, the 

Weymann et al. 1991). In the traditional dehnition of Wey- ,., ,. tttiat 111 1 

, , . , , optical contmua or LoBALs are more reddened compared 

mann et al. (1991), absorption troughs must be at least , tt-t-, a r .■ , 1 , ^ r. 

-1^ 1 J- c ; Qonn 1 -1 • to HiBALs, suggesting stronger dust extinction (e.g, Spray- 

2000 km s wide excluding the farst 3000 km s region , c ^^ nnn n ■ 1 j ^ 1 nnnn\ t v t -n a t 

, . , , , r berry & Foltz 1992; Reichard et al. 2003). In X-rays, LoBALs 

blue- ward from the emission lines to classify quasars as , , , . , , ,. , 1 ... tt- 

■niT^ri^T ■ ,n.. , n , ,r ^'^^o havc highcr gas absorption column densities than Hi- 

BALQbOs. Less strict dehnitions have also been used, tor -n a t / ^ 1 nnn-, ri n 1 ^ 1 mi 

, , r , , , BALs (e.g.. Green et al. 2001; Gallagher et al. 2002). There- 
example with a requirement or a trough to be at least „ ttiat 1 ^ i^- 

f 1 , / ^ , X tore, LoBALs probe a relatively obscured quasar population. 

1000 km s"^ wide (e.g., Trump et al. 2006). BALQSOs can ' ^ ^ h f p 

also be further divided into a population containing absorp- 



The origin of a small LoBAL fraction in quasars is unclear, 
, . , and it has been attributed to geometric effects (e.g., Elvis 

tion troughs from only high-ionization state species (e.g ^^^^^ evolutionary effects (e.g., Voit et al. 1993), 



C IV and N v; HiBALs) and a population that exhibits ab- 
sorption troughs in low-ionization species (e.g., Mg 11 and 



like the BALQSO population in general. There are several 
tentative arguments supporting the view that LoBALs are 
young quasars at a stage of blowing out obscuring materials. 
E-mail- dai@nhn ou edu First, several early studies of LoBAL fractions in the infra- 
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red band showed larger fractions (e.g., Boroson & Mayers 
1992) and associations with ultra- luminous infra-red galax- 
ies (e.g., Lipari et al. 1994; Canalizo & Stockton 2000). Sec- 
ond, a few optical spectral analyses suggested that the cov- 
ering fraction of the LoBAL wind is large (e.g., Voit et al. 
1993; Casebeer et al. 2008). Third, some radio spectra of 
LoBALs resemble those of compact steep spectrum or giga- 
hertz peaked spectrum sources, which are also candidates for 
young quasars (e.g., Montenegro- Montes et al. 2008; Liu et 
al. 2008). In particular, the LoBALs containing Fe absorp- 
tion troughs (FeLoBAL) are viewed as the most promising 
candidates for young quasars (e.g., Lipari et al. 2009). The 
fractions of LoBALs are important constraints on the origin 
of the LoBAL populations. 

Before studying the intrinsic fractions of LoBALs, it 
is important to compare the measurements of the intrinsic 
fractions of BALQSOs in the quasar population. Recently, 
a series of studies emerged on this topic. Dai, Shankar, 
and Sivakoff (2008a) studied the SDSS BALQSO (Trump 
et al. 2006) fractions in the 2MASS bands (Skrutskie et 
al. 2006), finding that the BALQSO fractions in the near- 
infrared (NIR) are twice those found in the optical band. In 
particular, we found the BALQSO fraction to be 20±2 per 
cent for the traditional BALQSOs that satisfy the stricter 
Weymann et al. (1991) definition and 43±2 per cent for the 
relaxed definition of Trump et al. (2006) that requires less 
broad absorption troughs. Dai et al. (2008a) argued that the 
BALQSO fractions measured in the NIR bands are closer to 
the intrinsic fractions, based on the observations that signifi- 
cant obscuration is associated with BALQSOs in the optical 
bands (e.g., Reichard et al. 2003), confirming the earlier esti- 
mate of Hewett & Foltz (2003). This result was confirmed by 
several studies, such as Ganguly & Brotherton (2008) using a 
different SDSS BALQSO catalogue, Maddox et al. (2008) us- 
ing the deeper UKIDSS survey, Shankar et al. (2008a) in the 
radio band, and Knigge et al. (2008) by correcting the frac- 
tion in the optical bands directly. In particular, Ganguly & 
Brotherton (2008) extended the study to include narrow and 
associated absorbers and found the overall outflowing AGNs 
to be 60 per cent of the total quasar population. Recently, 
Allen et al. (2011) claimed an even larger intrinsic fraction 
for the traditional BALQSOs of 41 per cent, by including 
the additional fraction of missing quasars that do not en- 
ter the SDSS survey. In particular, Allen et al. (2011) found 
that the completeness for BALQSOs and non-BALQSOs in 
SDSS is very similar at z < 2.1 and z > 3.6, but can differ 
at other redshifts, e.g., for z ~ 2.6 and z ~ 3.5. 

The larger fraction of BALQSOs makes the AGN wind 
a more promising candidate responsible for the feedback en- 
ergy that is needed to explain the co-evolution between black 
holes and their host galaxies. Understanding evolutionary 
versus geometric models of AGNs can not only probe the 
AGN feedback, but also constrain the nature of the feedback, 
whether it is kinetic from winds (e.g., Granato et al. 2004, 
Shankar et al. 2006, 2008b) or thermal (e.g., Di Matteo et 
al. 2005; Hopkins et al. 2006). If BALQSOs provide the ma- 
jority of the feedback energy, the feedback mechanism will 
be kinetic. Motivated by the results from BALQSOs, and 
the larger obscuration of LoBALs compared to HiBALs, we 
expect that the optical fractions for LoBALs are also biased 
low. This effect has already been noticed when only a few 
LoBALs were observed (Sprayberry & Foltz 1992); however, 



their study was limited by their small sample size. The large 
sample size enabled by SDSS warrants a new study on the 
intrinsic fractions of LoBALs. 

The radio properties of BALQSOs provide additional 
constraints on the nature of these objects. In an early study, 
Stocke et al. (1992) found no radio-loud BALQSOs within 
68 BALQSOs. Later studies showed that radio emission is 
present in BALQSOs (Francis et al. 1993; Brotherton et 
al. 1998; Becker et al. 2000); however, these BALQSOs are 
mostly radio-moderate. Matching the SDSS BALQSO cat- 
alogue in the FIRST survey (Becker et al. 1995), Shankar 
et al. (2008a) quantified the dependence of the BALQSO 
fraction on radio luminosities. We found that the BALQSO 
fraction drops at high radio luminosities confirming earlier 
claims of such an effect (Stocke et al. 1992; Becker et al. 
2001). In addition, Shankar et al. (2008a) also found that 
the BALQSO fraction at the low radio luminosity range is 
consistent with the NIR fraction of BALQSOs of Dai et al. 
(2008a). This result further supports the view that the NIR 
BALQSO fraction is close to the intrinsic fraction (mod- 
ulo corrections to the parent SDSS quasar selection), since 
there is also little absorption in the radio bands. The drop 
of the BALQSO fractions at high radio luminosities can be 
naturally explained under a geometric model of BALQSOs. 
If the radio emission has a preferred orientation, which is 
usually considered in the polar direction, the drop indicates 
that BALQSOs are less frequent in these viewing angles. Us- 
ing a unification model between radio-loud and radio-quiet 
quasars (e.g., Urry & Padovani 1995), we were able to suc- 
cessfully reproduce the trend, thus explaining the major- 
ity of the BALQSOs with radio emission under a geometric 
model. Exceptions still exist, such as the polar radio-loud 
BALQSOs, which were identified based on the radio vari- 
ability that implies too large brightness temperatures un- 
less the radio emission is relativistic (e.g., Zhou et al. 2006; 
Ghosh & Punsly 2007). Cold polar outflows are also present 
in blazars (e.g., Dai et al. 2008b), and they could be related 
to polar BALQSOs as the outflow continuously extends close 
to the polar axis. However, these objects are rare, and their 
implication for the total BALQSO population is still uncer- 
tain. 

In this paper, we study the intrinsic fraction of LoB- 
ALs by correlating SDSS quasars with detections in the NIR 
and radio bands. We also explore their radio properties and 
compare with BALQSOs to test whether LoBALs can be 
explained under a geometric or evolutionary model. We as- 
sume that 7^0 = 70 km s"^ Mpc"\ fim = 0.3, and = 0.7 
throughout the paper. 



2 SAMPLE SELECTION 

We select quasars from the SDSS DR3 quasar catalogue 
(Schneider et al. 2005) and the SDSS-BALQSO catalogue 
(Trump et al. 2006, T06 hereafter), where the authors also 
separately tagged the LoBAL population. In particular, we 
study quasars in the redshift range of 0.5 ^ 2: ^ 2.15, where 
the Mg 11 absorption line falls in the observed-frame optical 
band pass. Note that, in this redshift range, the complete- 
ness of BALQSOs and non-BALQSOs is very similar for the 
SDSS survey (Allen et al. 2011). For the FeLoBAL sample, 
we study the quasars in the redshift range of 1.19 ^ z ^ 2.24, 



© 2010 RAS, MNRAS OOO.IlVITOl 



Fractions of Low- Ionization BALQSOs 3 



13 



CP 14 
o 

}^ 15 F 

{/) 
(/) 
< 

^ 16 



17 
21 





' ' ' 1 1 


' ' 1 1 1 

■ A .-A'* ■ 


r 


A 
A 














Non-LoBAL -. 




■ ■ ■ ' I 


■■■ Bi-LoBAL A : 
AiNB-LoBAL n : 



20 19 18 17 16 
SDSS i (mog) 



15 




16 15 14 

2MASS Kg (mog) 



£ 



in 
c 

(D 
Q 

X 
3 



X 

o 



10" 
10^ 

10^ 

10' 

10° 



- 10-' 





500 : 


mog) 


400 'r 






d 
II 


300 r 


<u 

N 


200 ^ 


ISUI 






100 r 


z 


n ^ 



■4 r-< 



21 20 19 18 17 
SDSS i (mog) 



16 



non-LoBAL 
Bi-LoBAL*30 
AiNB-LoBal*30 



I : I 
I : i 





2 3 4 5 

SDSS r - 2MASS Kg (mog) 



15 



Figure 1. (top left) 2MASS Ks magnitude versus SDSS i magnitude for SDSS DR3 QSOs that are detected in all of the J, H, and 
Ks bands in the redshift range of 0.5 2: 2.15. The subsamples of QSOs that are not LoBALs (non-LoBALs; black dots), QSOs that 
satisfied the traditional Weymann et al. definition (BI- LoBALs; red triangles), and QSOs that satisfied the relaxed BAL definition of 
Trump et al. (2006), but did not satisfy the traditional definition (AINB-LoBALs; green squares) are displayed separately, (top right) 
SDSS i magnitude versus 1.4 GHz flux density, (bottom left) 2MASS Ks magnitude versus 1.4 GHz flux density, (bottom right) SDSS r— 
2MASS Kb colour for the three samples, where the histograms for the two LoBAL samples are multiplied arbitrarily by 30 for clarity. The 
LoBALs are significantly redder than the non-LoBAL population. The inset in the bottom right panel shows the cumulative distributions 
of the colour for the three samples. The K— S test results indicate that both of the LoBAL samples differ from the non-LoBAL sample 
with significances greater than 99.996 per cent, and that the two LoBAL samples are not significantly different from each other. 



which is the redshift range of FeLoBALs identified in TOG. 
Following Dai et al. (2008a) and Sliankar et al. (2008a) , we 
matcli the quasars using a 2" radius to entries from the full 
2MASS release, which extends to fluxes below the ofBcial 
completeness level of the 2MASS All-Sky Point Source Cat- 
alogue. To ensure that the matched database entries repre- 
sent detections, we require an in-band detection (rdJlg != 
0) that there is no confusion, contamination (cc_flg —— 0), 
or blending (bLflg ^ 1) in the source, and that no source is 
near an extended galaxy (gaLcontam —— and extjtey is 
null). The 99 per cent completeness levels of the database 
are J = 16.1, H — 15.5, and Ks = 15.1 mag. Recently, 
several new BALQSO catalogues have been published using 
more recent SDSS data releases (Gibson et al. 2009; Scaringi 
et al. 2009; Allen et al. 2011). In particular, the catalogue 
of Gibson et al. (2009, G09 hereafter) also provides identifi- 
cations of traditional Mg 11 LoBALs in the redshift range of 
0.55 ^ 2 ^ 2.15. We also briefly present the analysis to this 



sample as a comparison. We mainly present results based on 
the TOG sample, unless mentioned otherwise. 

Following Shankar et al. (2008a, 2010), we build a fuU 
FIRST-SDSS cross-correlation catalogue, containing all the 
detected radio components within 30"of an optical quasar. 
As our reference, following Schneider et al. (2005), we pri- 
marily present results for the FIRST-SDSS catalogue with 
radio counterparts identified within 2". We also test our re- 
sults by enlarging the radio matches to 5", and the results 
are consistent with those from the 2" matching. In addi- 
tion, many of the optical sources in SDSS are associated 
with more than one radio component in FIRST within 30", 
as expected if these sources are extended with jets and/or 
lobes separated from the central source. We use the sum of 
the flux densities as a proxy for the total radio luminosity 
of these sources. We require that at least one component is 
within our matching radius for the cases of multiple source 
matches. 

In the rest of the paper, we refer to LoBALs that sat- 
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Figure 2. Same plots as Fig. [T] but for FeLoBals in the rcdshift range of 1.19 ^ ^ ^5 2.24. The FeLoBALs are even more significantly 
redder than the non-LoBAL population. The K-S test result indicates that the r — Ks colours of the FcLoBAL population differ from 
those of the non-LoBAL population with a significance greater than 1 — 10~*. 



isfy the traditional definition of Weymann et al. (1991) as 
"BI-LoBALs", to LoBALs that are selected from TOG us- 
ing the relaxed AI definition as "AI-LoBALs", to LoBALs 
that satisfy the relaxed AI definition but not the traditional 
definition as "AINB- LoBALs" , and to the rest of the sam- 
ple without LoBAL features as "non-LoBALs" . Some studies 
have found that the AINB sample has properties besides the 
absorption troughs which differ from the BI sample, which 
may indicate a separate population (e.g., Knigge et al, 2008, 
Shankar et al. 2008a). We note that at the redshift ranges 
we study, the SDSS spectra cannot separate all the HiBALs 
from non-BALs because the high-ionization troughs are not 
redshifted into the optical spectral range. We also refer to 
the "2MASS sample" as quasars detected in all of the J, H, 
and Ks bands, and to the '^Ks complete sample" as quasars 
with Ks < 15.1 mag. We test our results using the two sam- 
ples, and they are not qualitatively different. Because of the 
scarcity of LoBALs, we choose to mainly present the results 
from the 2MASS sample to reduce the uncertainties of our 
fraction measurements. 



3 RESULTS 

3.1 Distributions of LoBALs in the Optical, NIR, 
and Radio Bands 

In three panels of Fig. [T] we plot the SDSS i magnitude, 
2MASS Ks magnitude, and FIRST 1.4 GHz fiux density 
against each other for three samples, BI-LoBALs, AINB- 
LoBALs, and non-LoBALs. The distributions of LoBAL 
samples are different from the non-LoBAL sample in all 
three cases. In the Ks versus i magnitude plot (Fig. [T] top 
left), the quasars are concentrated in a linear relation be- 
tween the i and Ks magnitudes with scatter. We can see 
that the LoBALs in the plot are systematically redder than 
the non-LoBAL population, as expected since the optical 
spectra of LoBALs show extra dust extinction compared 
to HiBALs and non-BALs (e.g., Reichard et al. 2003). We 
note that some of the scatter is due to quasar variability, 
since the SDSS and 2MASS data are taken from different 
epochs. However, this should not affect the mean colour dif- 
ference between the LoBALs and non-LoBALs in the plot. 
In the two panels involving radio flux in Fig. [T] we clearly 
see that the LoBALs are less populated in the high radio 
fiux regime. In the bottom-right panel of Fig. [T] we show 
the histograms of the SDSS r-2MASS Ks colour for the 
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BI-LoBALs, AINB-LoBALs, and non-LoBALs in our sam- 
ple as another example to demonstrate the redder colour of 
the LoBAL population. We also show the cumulative dis- 
tribution of the r — Ks colour in the inset, and perform 
the Kolmogorov-Smirnov (K-S) test on the distributions. 
The K-S test results show that the AI-LoBALs of our sam- 
ple are significantly different from the non-LoBALs with a 
very small probability (3 x 10~^^) for the two populations 
to arise from the same parent distribution. The situation 
is the same when we compare the sub-populations of the 
AI-LoBALs (BI-LoBALs and AINB-LoBALs) with the non- 
LoBALs, where we obtain K-S probabilities of 4 x 10~® and 
4 X 10"*. Comparing the BI-LoBALs and AINB-LoBALs 
in the plot, we find that the two colour distributions are 
not significantly different with a K-S probability of 0.28. In 
Fig.[2l we compare the properties of the FeLoBALs and non- 
LoB ALs of our sample using the same axes as in Fig. [T] We 
conclude that FeLoBALs are even more reddened than LoB- 
ALs, as expected. The radio powers of the FeLoBALs in the 
plot are similar to those of the BI-LoBALs. Ideally, these 
tests should be performed in narrower redshift bins, since 
there is quasar colour evolution with redshift; however, we 
are limited by the small sample size, and defer this analysis 
to future studies. 

3.2 The Intrinsic Fraction of LoBALs 

We present the fraction of LoBALs (BI-LoBALs, AI- 
LoBALs, and FeLoBALs) measured in the SDSS and 2MASS 
g,r,i, z, J, H , and Ks bands in Figs. |3] and [S] (left). We 
also list the numbers and fractions of LoBALs in Table \T[ 
where we use Gehrels's statistics to estimate the uncertain- 
ties (Gehrels 1986). We first calculate the absolute mag- 
nitudes of the quasars in the SDSS and 2MASS bands 
correcting for the Galactic extinction, but not correcting 
for the obscuration in the BALQSOs. We then calculate 
the fraction in the Ks band with an absolute magnitude 
limit of ~ —29.5 mag. For the rest of the bands, 

we set the limits so that the differences between the band 
limits are the same as their mean colour differences, i.e., 
Mx,um - MY,iim = Mx - My. We find that the LoBAL 
fractions are increasing from blue to red bands, similar to 
the results for BALQSOs (Dai et al. 2008a). This trend can 
be naturally explained considering the significant obscura- 
tion in LoBALs compared to HiBALs and non-BALs (e.g., 
Sprayberry & Foltz 1992), and the implication is that the 
2MASS fractions are closer to the intrinsic fraction of these 
objects. Comparing the BI-LoBAL samples from T06 and 
G09, we find that overall trend is consistent between the 
two samples. However, there is a systematic offset between 
the two samples, where the T06 sample has higher fractions. 

Following Dai et al. (2008a), we test the obscura- 
tion model by combining the quasar luminosity function 
(Richards et al. 2005) and the spectral differences be- 
tween LoBALs, HiBALs, and non-BALs (Reichard et al. 
2003), Model I hereafter. We assume that LoBALs, HiBALs, 
and non-BALs have the same intrinsic luminosity functions 
(Richards et al. 2005) with the only difference being in 
their normalizations. Since LoBALs and HiBALs are ob- 
scured, their observed luminosity functions are horizontally 
shifted to be less luminous, where the shift corresponds to 
the mean obscuration caused by dust extinction and ab- 



sorption. We calculate the obscuration based on the SDSS 
composite spectral models (Reichard et al. 2003). The shift 
is smaller in the red bands and larger in the blue bands. Fi- 
nally, we calculate the model observed fraction of LoBALs 
in the SDSS and 2MASS bands, and compare with the ob- 
servations. The only free parameter is the intrinsic fraction 
of LoBALs. Please refer to Dai et al. (2008a) for the detailed 
model. This assumes that the differential selection function 
for LoBALs and non-LoBALs in SDSS only depends on the 
different colours due to extinction and absorption troughs, 
which we believe is reasonable given our selected redshift 
range and the results of Allen et al. (2011). The model re- 
sults are presented in Fig. [3] (left), where the dashed lines 
are the intrinsic fractions of LoBALs and the solid lines con- 
nect the model results for the observed fraction of LoBALs. 
The square and filled circle symbols are the observed AI- 
LoBAL and BI-LoBAL fractions from T06. We also show 
the observed BI-LoBAL fractions from G09 as a compari- 
son. The model produces acceptable fits to the data with 
X'/dof = 1.5 and 2.2 for BLLoBALs and AI-LoBALs, and 
measures intrinsic fractions of LoBALs of 4.0 ± 0.6 per cent 
for BLLoBALs and 6.4 ± 0.8 per cent for ALLoBALs. We 
note that the model fits the data from the bluer bands bet- 
ter than the redder bands, which is possible that there is 
a separate IR luminous LoBAL population contributing to 
the LoBAL fractions, which we will discuss in detail later. 
The FeLoBAL fractions also show an increasing trend from 
the bluer g band to the redder Ks band (Fig. [5] left). We 
list the model results for using a pure geometric model in 
Table [1 

We plot the ALLoBAL, BLLoBAL (G09) and FeLoBAL 
fractions as a function of NIR and optical luminosities in 
Figs. |3]and[S] (right). In general, we find that the fractions 
are larger in the redder bands than bluer bands, consistent 
with the results found above. We also find that there is an 
increasing trend in LoBAL fractions as a function of lumi- 
nosity, especially in the redder 2MASS bands. The bluer 
g band does not show such a trend. This trend, especially 
in the redder bands, is not observed in BALQSOs in gen- 
eral (Dai et al. 2008a). Besides the bluer bands, which are 
more affected by the obscuration, the other bands have sim- 
ilar slopes for the increase of LoBAL fractions as a function 
of luminosity. However, the data at the luminous end have 
large uncertainties. Early IR surveys of LoBALs (with sam- 
ple sizes ~10) typically found LoBAL fractions of > 10 per 
cent (e.g., Boroson & Meyers 1992), higher than our intrin- 
sic fraction. It is possible that these early IR surveys have 
higher flux limits that can only probe the most IR luminous 
population, where we also find a large LoBAL fraction of 
> 10 per cent. 

It is possible that the increasing fraction with lumi- 
nosity is a result of increasing signal-to-noise ratio (S/N) 
spectra for quasars, such that it is more likely to identify 
sources as BALQSOs. This effect has been discussed previ- 
ously (Knigge et al. 2008; G09; Allen et al. 2011). However, 
our sample only probes the most luminous quasars after ap- 
plying the 2MASS flux limits, and in this regime, the effect 
due to S/N is small. For example, in Figure 11 of Allen et al. 
(2011), the C IV BAL fraction is a constant for the most lu- 
minous quasars. Moreover, if the trend of increasing fraction 
with luminosity is really due to the S/N of quasar spectra, 
it indicates that the intrinsic LoBAL fractions are close to 
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Figure 3. (left) Observed LoBAL fractions in optical and infrared bands. The filled circles and open squares are the observed BI-LoBAL 
and AI-LoBAL fractions from T06, and the open circles are from G09, where we shift the G09 data slight to the right for clarity. The BAL 
fractions increase from the blue to red bands, and the fractions can be modelled assuming there is significant obscuration in LoBALs. The 
heavy dashed line is the model intrinsic fraction for BI-LoBALs of T06, and the light dashed line is for that of AI-LoBALs of T06. The 
solid lines are the modeling results for the observed fractions, (right) We add an additional component of LoBALs at high luminosities 
to our model, where we obtain better fits to the data. The dashed lines show the model intrinsic fraction for the geometric component 
of LoBALs. 
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Figure 4. (left) The observed AI-LoBAL fractions as a function of optical and infrared luminosities. The g,i,z, and J band data are 
shifted right horizontally according to the mean colour differences with respect to the Ks band. The fractions increase from the blue 
to red bands because of the selection biases, especially in shorter wavelength bands. We also find that the AI-LoBAL fractions increase 
from low to high luminosities, especially in the redder bands. The solid line shows the model intrinsic fractions from Model II, where the 
dashed and dotted lines are the geometric and evolutionary component, respectively, (right) Same plot but for BI-LoBAL fractions from 
the G09 sample. 



~18 percent. This is not consistent with our Model I pre- 
diction after taking the obscuration eflFects into account. In 
the rest of the paper, we assume that the S/N has a small 
effect on the observed fractions in the luminosity range of 
our sample. 

We add another LoBAL component to our geometric 
model (Model II), based on the larger observed LoBAL frac- 
tions at high IR luminosities and the large values from 
a pure geometric fit in Fig. [3] (left). We call this compo- 
nent the evolutionary component. We model the evolution- 
ary component as having a power-law luminosity function 
only at high luminosities, ^evo = ^o,evoL" only when 
Mks ^ —29.5. We fit the observed LoBAL fractions from 
the g to Ks bands using Model II, and obtain consider- 



ably better fits. We show the Model II fits in Figs. [3] (right) 
and [5] (left), and list the fitting results in Table [2l The in- 
trinsic fractions from the geometric component are 3.3±0.4, 
5.4±0.5, and L3±0.3 per cent for BLLoBALs, AI-LoBALs, 
and FeLoBALs, respectively, and they are all smaller than 
the intrinsic fractions obtained from a pure geometric model 
fit (Model I). The fractions for the evolutionary component 
are functions of luminosities, and the total intrinsic fractions 
are the sum of the two components. The parameters of the 
evolutionary component, a and "3?o,bvo, are unconstrained 
from our data. We list one set of possible solutions in Ta- 
ble [2] as an example. After obtaining the best fit models, we 
plot the model intrinsic fractions of LoBALs as a function of 
Ka band luminosity in Figs.|4]and[5](right), for AI-LoBALs, 



© 2010 RAS, MNRAS OOO.fTl-ITol 



Fractions of Low- Ionization BALQSOs 7 



< 

m 
o 




to 
o 



c 
o 



o 
o 



< 
m 
o 



0.100 



0.010 



□ g-2.6mog 
i-2.3mag 
A z-2.3mog 
O J- 1 .2mog 

O K 



0.001 




-27 -28 -29 -30 -31 
Absolute Kj mog 



-32 -33 



Figure 5. (left) Observed FeLoBAL fractions from the g to Ks bands. The triangles are the data. The dotted line is the best fit from 
a pure geometric model, and dot-dashed line shows the intrinsic fraction from this model. The solid line is the best fit from Model II 
(geometric plus evolutionary), and the dashed line shows intrinsic fraction for the geometric component, (right) Observed FeLoBAL 
fractions as a function of optical and infrared luminosities. The solid line is the model intrinsic FeLoBAL fraction from Model II, and 
the dashed and dotted lines are the geometric and evolutionary component, respectively. 



Table 1. Fractions of LoBALs in the 2MASS and SDSS Bands 





Sample Redshift Range 


u 


9 


r 


i 


z 


J 


H 


Ks 


Limiting Mag 




-26.7 


-26.9 


-27.1 


-27.2 


-27.2 


-28.3 


-28.9 


-29.5 



N (Total) 


T06 


0.50-2.15 


1460 


1181 


1062 


1312 


1518 


776 


724 


533 


N (BI-LoBAL) 






5 


11 


17 


19 


24 


26 


29 


22 


f (BI-LoBAL) 








9+°-'' 


'-■"-OA 


'-■"-0.3 


'-■"-0.3 


3 4+0.8 


^■'-'-0.7 


4 1+1-1 


N (AI-LoBAL) 






12 


19 


26 


30 


39 


38 


45 


36 


f (AI-LoBAL) 






o.8t«:i 




9 t;+0.6 


9 O+0.5 


'^■"-0.4 


4 9+°-^ 


6 2+1-1 


6.7tl-l 



N (Total) 1.19-2.24 1584 1366 1207 1486 1783 843 798 610 

N (FeLoBAL) 1 2 8 9 15 19 22 19 

f (FeLoBAL) O^OStZl 0-2to.? ^■'^-o.l 0-6to.2 O^^-oi ^■^t°o.l S.S^;^ S.l^-? 

N (Total) G09 0.55-2.15 2225 1817 1621 2030 2337 1252 1313 991 

N (BLLoBAL) 2 5 15 19 27 37 43 32 

f(BI-LoBAL) Om+_l:ll 0.3tri 0-9tg:l 0-91^2 l-^lg:^ 3.0lgj 3.3tg-j 3.2tg;^ 

Notes: The fractions arc in per cents. 

Table 2. Modeling Results for the Intrinsic Fractions of LoBALs Using the 2MASS and SDSS Data 





Model 


AI TOO 


BI TOO 


BI G09 


Fe TOO 


Intrinsic Fraction (per cent) 


I 


6.4±0.8 


4.0±0.6 


2.8±0.4 


1.9±0.4 


xVdof 




13.4/6 


9.0/6 


21.6/6 


16.9/6 


Intrinsic Geometric Fraction (per cent) 


II 


5.4±0.5 


3.3±0.4 


2.0±0.3 


1.3±0.3 


a 
















^0,EVo/^0,Total 




6.1x10-5 


4.2x10-5 


5.0x10-5 


4.9x10-5 


xVdof 




2.7/4 


1.5/4 


2.6/4 


2.4/4 



Notes: Model I is a pure geometric model, where the only free parameter is the intrinsic fractions of LoBALs. Model II is a two 
components model, where we add an additional component (evolutionary) to Model I. The parameter a is the power-law slope of the 
luminosity function of the evolutionary component for LoBALs, ^evO ^0 EVO^"'- The parameter ^o,EVO /^O.Total is the 
normalization ratio between the luminosity function of the evolutionary component for LoBALs and total quasars. Both a and 
^0,EVO /^0,Total '^^^ unconstrained from the data. We list one set of possible solutions that produce good fits in this table. 
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BI-LoBALs (G09), and FeLoBALs, respectively. In all three 
cases, we find that the models fit the observed Ks band frac- 
tions well. The models (solid lines), which are the intrinsic 
fractions, represent an envelope for observed fractions. The 
observed fractions are below the models since they are bi- 
ased with different degrees. The Ks band fractions should 
be closest to the model, since there is less absorption in 
the Ks band. The observed Ks band fractions for LoBALs 
fainter than Ks > —29 are slightly below the models, which 
is possibly associated with the incomplete Ks band data. 

3.3 Radio Properties of LoBALs 

We show the fraction of LoBALs as a function of the 1.4 GHz 
radio luminosity in Fig. |6] (left). In the 0.5 ^ 2 ^ 2.15 red- 
shift range, we find 2650, 95, and 49 matches in FIRST 
for the SDSS quasars, AI-LoBALs, and BI-LoBALs, respec- 
tively. In the 1.19 ^ ^ ^ 2.24 redshift range, we find 1621 
and 19 matches in FIRST for the SDSS quasars and FeLoB- 
ALs, respectively. We find that the fraction of LoBALs de- 
creases with increasing radio luminosities. This is true in 
all of the sub-samples of LoBALs (BI-LoBALs, AI-LoBALs, 
AINB- LoBALs, and FeLoBALs). We do not show the AINB- 
LoBAL fractions in Fig. |4] for clarity, but the numbers can 
be obtained by subtracting the BI-LoBAL fractions from 
the AI-LoBAL fractions. In addition, the LoBAL fraction 
at the low radio luminosity end is consistent with the intrin- 
sic fraction obtained in Section [3.21 This is not unexpected 
since there is also little obscuration in the radio band. These 
two features are similar to BALQSOs (Shankar et al. 2008a). 

The decrease in the LoBAL fractions as a function of 
radio luminosity can be easily explained under a geometric 
model for LoBALs, if we assume that LoBALs are viewed 
in the lines of sight close to the accretion disc, while the 
boosted radio emission is viewed close to the polar direc- 
tion. We test this model following Shankar et al. (2008a). We 
assume that the radio emission is composed of one weaker 
uniform component and another strong, beamed component 
at the polar direction. LoBALs are modelled to be within a 
solid angle close to the equatorial plane, where the open an- 
gle is a free parameter. The rest of the parameters, such as 
the slope of the luminosity function, Lorentz factor, maxi- 
mum beaming angle, are all fixed according to the literature 
(Urry & Padovani 1991; Urry & Padovani 1995; De Zotti et 
al. 2005; Richards et al. 2006; Jiang et al. 2007; Padovani 
et al. 2007). Please refer to Shankar et al. (2008a) for the 
detailed model. Using this model, we find good fits to the 
data with x^/dof = 1.4, 0.12, and 1.0 for BLLoBALs, AL 
LoBALs, and FeLoBALs, respectively, and the correspond- 
ing intrinsic fractions of these objects are 4.0±0.7, 8.4±1.0, 
and 3.6 ± 1.0 per cent in quasars. These fractions are con- 
sistent with the LoBAL fractions obtained using a pure ge- 
ometric model (Model I) in Section 3.2; however, they are 
higher than those obtained using a two components model 
(Model II). 

Since the infrared and optical data favour an additional 
evolutionary component for LoBALs, we try to incorporate 
this component in our radio model. From the radio flux ver- 
sus the Ks magnitude plots in Figs. [1] and (2] it is likely that 
the Ks bright LoBALs have relatively faint radio fluxes. 
Since there is little absorption in the radio band and our 
model shows that the measured fractions at low radio lumi- 



nosities are close to the intrinsic fractions, we can subtract 
the evolutionary component assuming the ratio of the evo- 
lutionary component to the geometric component is a con- 
stant at low radio luminosities. Based on this, we subtract 
the evolutionary component in the LoBAL fractions at low 
radio luminosities {Liaghz ^ lO^^ergs s"^ Hz"^) using the 
ratio between the model intrinsic fraction for the geometric 
component in Model II and the model intrinsic fraction in 
Model I in section Section 3.2. The resulting data are shown 
in Fig. |6] (right), where we flnd lower fractions compared to 
the left panel at Liaghz ^ lO'^'^ergs s~^ Hz~^. We flt these 
data using a pure geometric model as above, and flnd intrin- 
sic fractions of 3.6±0.6, 7.2±0.9, and 2.8 ± 0.8 per cent for 
the geometric components for BI-LoBALs, AI-LoBALs, and 
FeLoBALs, respectively. We also sightly improve the fltting 
statistics. The fltting results are listed in Tabled 



4 SUMMARY AND DISCUSSION 

We flnd signiflcantly high fractions of LoBALs in the quasar 
population compared to the values obtained using optical 
data only. For example, the BI-LoBAL and AI-LoBAL frac- 
tions in the optical data were measured as 0.55 and 1.31 
per cent (T06), while our results are 5-7 times larger. Al- 
though the final intrinsic fractions depend on the choice of 
catalogues, the overall trend is found to be the same. For ex- 
ample, we perform a similar analysis to the BI-LoBALs from 
the G09 sample, and also flnd large intrinsic LoBAL frac- 
tions. Although there is a systematic offset between the BI- 
LoBAL fractions from the T06 and G09 samples, the overall 
trend for the observed fractions is the same. Our intrinsic 
fractions are obtained using two independent methods, one 
from the NIR and optical data and the other from the radio 
data, and the results are mutually consistent between the 
two methods. Combining the estimates from the two meth- 
ods using the least square (minimum variance) method, we 
flnd that the intrinsic fractions for BI-LoBALs, AI-LoBALs, 
and FeLoBALs are 4.0 ±0.5, 7.1 ±0.6, and 2.1 ±0.3 per cent, 
respectively, using a pure geometric model. For our hybrid 
model with both the geometric and evolutionary compo- 
nents, which fit the data better, the corresponding intrinsic 
fractions for the geometric component are 3.4±0.3, 5.8±0.4, 
and 1.5 ± 0.3 per cent. The intrinsic fractions for the evolu- 
tionary component for LoBALs are functions of luminosities, 
and the total intrinsic fractions of LoBALs are the sums of 
the two components. The final combined estimates of the in- 
trinsic LoBAL fractions are listed in Tabled] The results are 
not unexpected considering the large obscuration observed 
in LoBALs (e.g., Sprayberry & Foltz 1992). Although we 
find significantly larger intrinsic fractions of LoBALs, they 
still represent a small portion of the total population. Com- 
pared with the intrinsic fractions of 20 ± 2 and 43 ± 2 per 
cent for BI-BALQSOs and AI-BALQSOs (Dai et al. 2008a), 
the LoBALs are about 20 per cent of BALQSOs. 

Our method of calculating the intrinsic fractions of 
BALQSOs still depends on the completeness of optical 
quasar surveys. The fraction of quasars that do not enter 
the optical surveys was not be accounted in this or our pre- 
vious papers. It was estimated that SDSS is about 90 per 
cent complete at 2; < 2.2 and i < 19.1 mag (Richards et 
al. 2002). The remaining 10 per cent of quasars, which are 
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Figure 6. (left) Fractions of LoBALs and FeLoBALs as a function of the 1.4 GHz radio luminosity. The data are slightly shifted 
horizontally for clarity. At the low luminosity end, the fractions are consistent with the intrinsic fractions modelled through the NIR 
and optical bands. The fractions are significantly smaller at high radio luminosities. Using a geometric model for LoBALs (solid, dotted, 
and dashed lines for BI-LoBALs, AI-LoBALs, and FeLoBALs, respectively), we are able to fit the trend and obtain intrinsic fractions of 
4.0±0.7, 8.4±1.0, and 3.6±1.0 per cent for BI-LoBALs, AI-LoBALs, and FeLoBALs, respectively, (right) We subtract the evolutionary 
component at Li ^quz ^ 10'^'^ ergs Hz~^, and fit the modified data using a geometric model. We find intrinsic fractions of 3.6±0.6, 
7.2±0.9, and 2.8 it 0.8 per cent for the geometric component for BI-LoBALs, AI-LoBALs, and FeLoBALs, respectively. 



Table 3. Modeling Results for the Intrinsic Fractions of LoBALs Using the FIRST Data 





Model 


ALLoBAL 


BI-LoBAL 


FcLoBAL 


Intrinsic Fraction (per cent) 


Model I 


8.4±1.0 


4.0±0.7 


3.6±1.0 


X'/dof 




0.9/7 


9.9/7 


6.3/6 


Intrinsic Geometric Fraction (per cent) 


Model II 


7.2±0.9 


3.6±0.6 


2.8±0.8 


X^/dof 




0.5/7 


7.5/7 


4.7/7 



thought to be highly obscured, could potentiaUy all be LoB- 
ALs or FeLoBALs. Thus, in most optimistic estimates, the 
fractions of LoBALs or FeLoBALs can reach ~15 per cent 
for the geometric component. However, the nature of the 
obscured quasars are still uncertain, and we are not sure 
whether the uhra-violet spectra of these quasars still show 
broad absorption lines, if the continua of these quasars are 
mostly obscured. Therefore, the fractions quoted in our pa- 
per represent conservative estimates based on observations. 

The LoBAL fractions in the radio band are particu- 
larly interesting, since we find that the LoBAL fractions 
decrease with increasing radio luminosities. This confirms 
the early result of Becker et al. (2000) with about a dozen 
LoBALs. The trend is similar to that found in the total 
BALQSO population (Shankar et al. 2008a). This trend 
found in both the total BALQSO population and LoBALs 
suggests that the majority of LoBALs and BALQSOs can 
be united under a similar physical scheme. In Shankar et 
al. (2008a), we favoured a geometric model to interpret 
the trend. Applying the geometric model of Shankar et al. 
(2008a) to LoBALs, we can successfully reproduce the data 
for BLLoBALs, ALLoBALs, and FeLoBALs. However, it 
is problematic under a pure evolutionary model to explain 
the radio- luminosity/LoBAL fraction trend as argued by 
Shankar et al. (2008a). The main question is why quasars 
spend the same fraction of time as BALQSOs in the radio- 
loud stage and radio-quiet stage. Therefore, we argue that 



the majority of BALQSOs and LoBALs can be understood 
in a geometric model. Polar BALQSOs/outfiows (Zhou et 
al. 2006; Ghosh & Punsly 2007; Dai et al. 2008b) present 
a challenge to our results; however, these objects are rare 
and we are uncertain about their implications to the total 
quasar population. A modification of the geometric model 
to have both disc and polar outflows (e.g., Borguet & Hut- 
semekers 2010) may be needed to incorporate these objects. 
We also note that LoBALs dominate this population of polar 
BALQSO candidates (Ghosh & Punsly 2007). 

There are other indications, such as the association with 
ULIRGs, radio spectra, and large covering fractions from 
spectral modeling, arguing that LoBALs belong to an ear- 
lier evolutionary stage of quasar population. However, most 
of these studies are based on a small sample size and may 
not extrapolate to the whole LoBAL population. Urrutia et 
al. (2009) studied the fraction of LoBALs in the dust red- 
dened quasars at high redshift, finding that all except one 
are LoBALs, supporting the young nature of LoBALs. How- 
ever, the authors also noted that their selection method may 
be biased favouring LoBALs since they are associated with 
large dust reddening. In our study, we find LoBALs and 
BALQSOs are similar in most aspects, except that the frac- 
tion of LoBALs increases with increasing NIR luminosities. 
This is not consistent with BALQSOs in general, because 
the BALQSO fractions are mostly constant with increasing 
NIR luminosities for ALBALQSOs (Dai et al. 2008a). At the 
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Table 4. Final Intrinsic Fractions of LoBALs by Combining the Two Methods 





Model 


AI-LoBAL 


BI-LoBAL 


FeLoBAL 


Intrinsic Fraction (per cent) 


Model I 


7.1±0.6 


4.0±0.5 


2.1±0.3 


Intrinsic Geometric Fraction (per cent) 


Model II 


5.8±0.4 


3.4±0.3 


1.5±0.3 



NIR luminous end, the observed LoBAL fraction is higher, 
although with large uncertainties, than the intrinsic fraction 
that we obtain for a pure geometric model. 

It is possible that a portion of NIR luminous LoBALs 
are special compared to the rest of the population, and at an 
early evolutionary stage of the quasar cycle. This will rec- 
oncile some observations supporting the young quasar inter- 
pretation for LoBALs. We add an evolutionary LoBAL com- 
ponent to our model, and generally obtain a better fit to the 
data. However, the parameters of this evolutionary compo- 
nent are unconstrained from the data. The combination of 
a short early evolution model, with a large covering fraction 
close to 100 per cent, and a subsequent, longer geometric 
model, with covering fractions consistent with the intrin- 
sic fraction measured in this paper and Dai et al. (2008a), 
could simultaneously account for many of the differences 
and similarities between LoBALs and HiBALs. The longer- 
lived geometric model mainly sets the intrinsic fractions of 
BALQSOs of various species and the measured fraction as 
a function of their radio luminosities. If there is some small 
spherical outflow component at early times, this might also 
explain the predominance of LoBALs among the rare po- 
lar BALQSOs. Deeper IR and radio surveys are needed to 
increase the sample size and confirm this claim. 
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